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The coincidence problem for three-dimensional discrete structures with icosahedral symme- 
try is reinvestigated. We present a parametric description of the coincidence rotations based 
on special quaternions, called icosian numbers. In particular, we give a characterization of the 
possible coincidence indices (E-factors) and present a complete list of the possible rotations 
with E < 100. 

1 Introduction 



^ ■ The symmetries of crystallographic lattices are well known: these are the operations that map the 

whole lattice onto itself, so that each lattice point coincides exactly onto another one. The symmetry 
operations are translations, rotations, reflections and combinations thereof. It is, however, also 
possible to map only a part of the lattice points onto the others: a sublattice emerges. If the 
sublattice has full rank (i.e. same dimension as the original lattice), but possibly lower symmetry, 
' it is called a coincidence site lattice, or CSL for short. The CSL is characterized by its index in the 

ON ■ original lattice. This number is also called the E-factor. CSLs are useful for several aspects. 

First of all, there are mathematical reasons: the point symmetry group of the lattice Q is a 

S subgroup of the symmetry group of the isotropic Euclidean space, SO(n, IR). The set of coincidence 
rotations also forms a group, and includes the group of symmetries. In particular, it contains 
. rotations with rational angle. Therefore we have SO(n,1.) C CSL — group C SO(n, IR), and it 

would be desirable to know more about the structure of the CSL group. Since the point symmetry 
group may be a maximal finite subgroup of SO(n, IR), the CSL group must have infinite order. 

Second there are crystallographic reasons: A real crystal consists of several grains. It has been 
shown that these components are mostly not oriented arbitrarily with respect to each other, but 
that two adjacent grains often have a common coincidence lattice with low index in the full lattice. 
The CSL are important in crystallography because they allow a nontrivial classification of grain 
boundaries and because small-unit-cell CSL grain boundaries seem to be energetically favoured []. 
It has been shown that coincidences of vertices appear also in quasicrystallinc tilings [] . Quasicrystals 
also have grain boundaries like ordinary crystals, and one should therefore know the coincidence 
site quasilattices []. Twinning has been observed in experiment [], as well as interfaces between 
icosahedral and decagonal quasicrystals []. Similar to the crystal the CSL grain boundary should 
be energetically favoured. To understand the occuring coincidence indices one has to enumerated 
the possible CSL. The grain boundary between crystals itself may be quasiperiodic as predicted by 
Rivier and Lawrence []. Indirect evidence was provided already in the 1970s [], and later observed 
in the growth of quasicrystalline grain at the grain boundary between two crystals [] . 

Some microstructures with icosahedral (pseudo) symmetry are composed of 'approximants' by 
twinning and multiple twinning. The grain boundaries of these nanodomains will also be low energy 
grain boundaries, but their orientation relation is different from those ones described in the present 
paper. 

The first one who has dealt with generating functions and the enumeration of CSL seems to have 
been Ranganathan [] . The classification of the crystallographic CSL has been obtained by several 
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authors in the '70's and '80's [], and there exist tables of cubic and hexagonal CSL. But there 
has been no systematic mathematical treatment of the problem until recently, see [] and references 
therein. 

The interest in CSL was renewed with the discovery of the quasicrystals. It has been found 
that this new type of materials also exhibits multiple grains, twin relationships and coincidence 
(quasi-)lattices. 

In a quasicrystal the number of generating basis vectors is larger than the space dimension. This 
means that there exists no lattice with a minimal distance in physical space but a dense module. 
The CSL has to be replaced by a coincidence site module (CSM). A quasilattice is generated if only 
certain points are selected from the module by a window function. 

Recently the complete classification of all possible CSL / CSM for crystals and quasicrystals with 
rotational symmetry of any order and an analysis of the structure of the CSL/CSM group has been 
achieved by Pleasants et al. [] with the help of number theoretical methods. This means that the 
CSL/CSM case is solved completely in two dimensions. In three and four dimensions, similar (but 
less complete) results are known through the use of quaternions, see Baake and Pleasants [] and 
Baake []. 

In three dimensions the problem has also been adressed by Warrington and Luck [], who enu- 
merated the icosahedral CSM with small index, and by Radulescu and Warrington [] . In this paper 
we will generalize the systematic and appealing treatment of the crystallographic CSL in three di- 
mensions by Grimmer [] to icosahedral quasicrystals. The method is suitable for an implementation 
in a computer program and allows an easy enumeration of all CSM up to high indices. Our main 
purpose is therefore the presentation of complete tables of icosahedral CSM representatives. 

2 Formulation of the problem 

A coincindence rotation C maps a lattice L onto a lattice L': 

L' = CL (1) 

The rotation C is not unique, since we can apply symmetry rotations R and S to the lattices L = RL 
and L' = SL' , resulting in a conjugation of C": 

SL' = CRL or L 1 = S~ l CRL (2) 

By comparison with Eq. |l| we may introduce a new coincidence rotation C: 

C = S- X CR (3) 

This means that we have to count the coincindence rotations C with respect to double coset classes: 

(S,R) £ SO (3) x SO(3) £ SO(4)/Z 2 

First of all we need a proper parametrization of the rotation matrices. A three-dimensional matrix 
has nine entries, but only three are independent: the axis (of length one) and the angle. Cayley's 
parametrization through quaternions is useful here []. If q = (k, A,/i, v) is the quaternion, the first 
entry n parametrizes the angle: 



cos(0) 



*-2 \2 ,,2 ,,2 

K — A — [I — V 



K 2 + A 2 + [l 2 + V 2 

and the vector (A,/i, v) represents the rotation axis. The rotation matrix is therefore given by: 

1 I k 2 + A 2 - [i 2 - v 2 -2kv + 2\fi 2k{i + 2\v 

R = — 2kv + 2Xfi k 2 - A 2 + [i 2 - v 2 -2kA + 2[iv 

\ -2k[i + 2\v 2k\ + 2{iv k 2 - A 2 - \x 2 + v 2 



Figure 1: The asymmetric unit of the four-dimensional polytope {5,3,3}. The dash-dotted line is 
the 5-fold axis, the long-dashed lines are 2-fold axis and the short-dashed lines 3-fold axis. 



a1 = xa2+T a3 a0 = T 2 a l + Ta2 




a2=0 



The norm of a quaternion is: 

|q| 2 = K 2 + A 2 +,i 2 + ^ 

We have a two-to-one homomorphism of the quaternions and the rotation matrices since i?(q) = 
R(— q). In the quaternionic formulation of the coincidence problem we get [] (small letters indicate 
quaternions): 

L = C'L -> l^c'lc'- 1 

C = S^CR -► c' = s^cr (4) 

Similar to the three-dimensional Cayley parametrization there exists also a four-dimensional one: 
Due to the homomorphisms between the quaternions, SU(2) and 50(3) and the homomorphism 
between 5*0(3) x 50(3), SU{2) x SU{2) and 50(4), we can parametrize four-dimensional rotation 
matrices by a pair of quaternions [] . The equation is thus further transformed into 

c' = s J cr C = MC (5) 

Now C represents a four-dimensional vector and M a four-dimensional rotation matrix parametrized 
by two quaternions r = (k, I, m, n, ) and s = (a, 6, c, d). The rotation matix has the form: 

( ak + bl + cm + dn —al + bk + cn — dm —am — bn + ck + dl —an + bm — cl + dk \ 

al — bk + cn — dm ak + bl — cm — dn —an + bm + cl — dk am + bn + ck + dl 

am — bn — ck + dl an + bm + cl + dk ak — bl + cm — dn —al — bk + cn + dm 

V an + bm — cl — dk —am + bn — ck + dl al + bk + cn + dm ak — bl — cm + dn J 

What is the meaning of this transformation? Equation [s] is nothing else than the symmetry 
description of a four-dimensional polytope! A detailed description of the polytope may be found in 
Ref. []. We have transformed the problem of enumerating double coset classes into the problem of 
analyzing a four-dimensional polytope. If we calculate the orbit of a certain point of the polytope 
we find that the size of the orbit gives the number of equivalent CSM, and the type of the orbit tells 
us the symmetry of the CSM. Practically, we have to find a standard representation for the CSM. 
This is the disorientation [] : within a class of equivalent coincidence rotations there is one where 
the rotation angle transforming L into L' is minimal. Therefore k in q is maximal. Translated into 
the four-dimensional language of the polytopes we get a set of linear inequalities [] which define an 
asymmetric unit. Fig. [I] displays the icosahedral case. The crystallographic cases can be found in 
Ref. []. 



In order to count the orbits we have to fill the polytope with points. In the cubic case it has 
been shown [] that it is sufficient that the entries of the quaternion q = (k, A, u, v) are integers. In 
the icosahedral case we have to use numbers of the type m + nr with m and n integers and r the 
golden mean + l)/2 []. These numbers are elements of the ring of icosians generated by the 
quaternions (1,0,0,0), 1/2(1,1,1,1) and 1/2(t, 1, — 1/r, 0). These quaternions generate the group 
Y of order 120, and the ring of icosians I are all integer linear combinations. 

There exists a further interesting module: It consists of all points in Z[r], has cubic symmetry 
and contains the icosahedral modules as subsets (for details see []). The coincidence rotations are 
generated by the group 5*0(3, Q(r)), and the quaternions have entries of type m + tit without any 
restriction. 

3 Technicalities 

To avoid counting equivalent quaternions several times we have to demand that three conditions are 
fulfilled: First, we take only quaternions whose components have no common factor: gcd(«, A, /i, v) — 
1. Second, the entries of q are ordered with decreasing size. Third, the norm of a number m + nr 
[]is: 

N(m + nr) = mi 2 + run — n 2 

This means that 

N(±r k ) = (-l) fc 

with k an integer. To take this last condition into account it is convenient to restrict the smallest 
non-zero quaternion component between 1 and r. 
The index for an icosahedral CSM is given by: 

Z = N(\q\ 2 )/2 e 

where i = if one of the (coprime) quaternion components is odd and t — 2 if all are even. For the 
cubic case it is: 

E = |iV(dcn(i?))| 
where den(i?) is the denominator of the rotation matrix. 

4 Results 

It is easy to count icosahedral CSM and the CSM of Z[r] on a workstation by looping through the 
unit icosians or through the quaternions with m + nr entires using the algorithm described above. 
The correctness of the results can be checked by the generating function for the number of CSM 
#(E) given in Ref. []. We have tabulated the CSM with a rotation axis parallel to a symmetry axis 
up to E = 5000, with an axis in a general up to E = 500. Some special cases of the order E = 20000 
have also been checked. Summary tables of the results are given in Tab. representatives for the 
icosahedral case in Tabs. HJ|,||. The cubic case is presented in Tabs. |l],|]|j, an d [7|,pp|JTo|. 

It is interesting to note that the number of CSM for the cubic and icosahedral case are identical 
if E is odd, but if E is even, iV cufc (E) = 8 * iV lco (E/4). 

Miscellany 

The results described in this paper have been presented at the XXI. Int. Coll. on Group Theoret. 
Methods in Physics, July 15th-19th, in Goslar, Germany. The tables listed in this paper may be 
obtained from the authors through email from johannes.roth@itap.physik.uni-stuttgart.de. 
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Tabic 1: Summary of the CSM orbits. I. The part left of the double line are the results for the 
icosahedral case, the part to the right is for the cubic module Z[t] 3 . Cn denotes a CSM axis parallel 
to n-fold axis, M a rotation axis in a mirror plane and G a rotation axis in a general position. The 
numbers #(cr) give the sum of the lengths of the orbits. Entries of type a*b denote a*b CSM with 
b algebraic conjugates. The symmetry groups for £ < 10 are explicitly given. Primes indicate a 
CSM occuring in a standard cubic lattice. 
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Table 2: Summary of the CSM orbits. II. Notation as in Tab. |. 
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Table 3: Summary of the CSM orbits. III. This table lists the coincidence rotations for the cubic 
module Z[r] 3 with even £. Notation as in Tab. [|. 
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Tabic 4: Representatives for coincidence rotations with S < 100 for the icosahcdral module. I. The 
two sections indicate rotation axis along 5- and 3-fold symmetry axis, respectively. The coordinates 
are given as m + rn. The first column lists the conjugate pairs and the last column the disorientation 
angle. 
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Table 5: Representatives for coincidence rotations with X < 100 for the icosahedral module. II. 
Same as Tab. for rotations along 2-fold symmetry axis. 
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Table 6: Representatives for coincidence rotations with £ < 100 for the icosahedral module. III. 
Same as Tab. for rotations in mirror planes or general positions (last entry only). 
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Tabic 7: Representatives for coincidence rotations with odd X < 100 for the cubic module Z[t] 3 . I. 
The coordinates are given as m + rn. The sixth column gives the orbit length. The last column 
lists the disorientation angle. 
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Table 8: Representatives for coincidence rotations 
Notation as in Tab. R. 
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< 100 for the cubic module Z[r]. II. 
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Table 9: Representatives for coincidence rotations with odd £ < 100 for the cubic module Z[r]. 
Notation as in Tab. 0. 
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Tabic 10: Representatives for coincidence rotations with even £ < 100 for the cubic module Z[r]. IV. 
Notation as in Tabl. 0. 
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